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Abstract. We show that a quantum control procedure on a two-level system including 
çSJ , dissipation gives rise to the semi-group corresponding to the Lie àlgebra gl(3, R) © i? 3 . The 

■ physical evolution may be modelled by the action of this semi-group on a 3-vector as it moves 

q ' inside the Bloch sphere, in the Bloch ball. 

1. Introduction 

> 

Recent developments in quantum computing have emphasized the need for a realistic analysis 
of dissipation in systems which have the potential for use as qubits. In this note we discuss 
the effects of control and dissipation on a two-level system. For a single qubit pure state, it 
is well known that the unitary evolution may be visualised as the movement of a vector, the 
Bloch vector, on the surface of a 2-sphere, the Bloch sphere. In this note we extend the idea 
to a two-level mixed state. For this system, unitary evolution is on a spherical shell within the 
Bloch Sphere. Dissipation causes more general motion within the Bloch ball. This motion 
corresponds to the action of a certain semi-group. 

We also show that the effects of dissipation may not be compensated by the interaction 
with the external control. However, taking into account the effects of measurement, which 
may be modelled by certain projection operators, the dissipative effects may indeed be 
modiíied, allowing more effective control of the system. 
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2. Dynamics of dissipative quantum control systems 



In pure-state quantum mechanics the state of the system is usually represented by a 
wavefunction \^>), which is an element of a Hilbert space H. For dissipative quantum systems, 
however, a quantum statistical mechanics formulation is necessary since dissipative effects 
can and do convert pure states into statistical ensembles and vice versa. In this case, the state 
of the system is represented by a density operator p, whose diagonal elements determine 
the populations of the energy eigenstates, while the off-diagonal elements determine the 
coherences between energy eigenstates, which distinguish coherent superposition states 
1^) = Hn=i c n\ n ) from statistical ensembles of energy eigenstates (i.e., mixed states) 
p = Y,n=i w n \n)(n\. For a non-dissipative system the time evolution of the density matrix 
p(t) with p(t ) = p is governed by 

pit) = ü(t) Po ü(ty, (d 

where U (í) is the time-evolution operator satisfying the Schrodinger equation 

ihjÜ(t) = H(ï)Ü(t), 17(0) =1, (2) 
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where J is the identity operator. p(t) also satisfies the quantum Liouville equation 

ihj t p(t) = [H(f),p(t)\ = È{ï)p{t) - p(t)H(í). (3) 
H{i ) is the total Hamiltonian of the system, which depends on a set of control fields f m ; 

M 

H(f)=È +Y,fm{t)H m , (4) 

m=l 

where H is the internal Hamiltonian and H m is the interaction Hamiltonian for the field f m 

for 1 < m < M. 

The advantage of the Liouville equation @ over the unitary evolution equation (HJ) is 
that it can easily be adapted to dissipative systems by adding a dissipation (super-)operator 

Mp(í)]: 

M 

ihp(t) = [H ,p(t)\ + £ f m {t)[H m ,p(t)) + iKL D \p(t)]. (5) 

m=l 

In general, uncontrollable interactions of the system with its environment lead to two 
types of dissipation: phase decoherence (dephasing) and population relaxation (decay). The 
former occurs when the interaction with the enviroment destroys the phase correlations 
between states, which leads to a decay of the off-diagonal elements of the density matrix: 

Pknit) = -~([È(f),p(t)]) kn - T kn p kn (t) (6) 

where T k n (for k ^ n) is the dephasing rate between \k) and \n). The latter happens, for 
instance, when a quantum particle in state \n) spontaneously emits a photon and decays to 
another quantum state \k), which changes the populations according to 

Pnn{t) = -~([H(f),p(t)]) nn + hnkPkkit) - ^knPnnif)} (V) 

k^n 

where ^knPnn is the population loss for level |n) due to transitions |n) — > \k), and jnkPkk 
is the population gain caused by transitions \k) — > \n). The population relaxation rate jkn 
is determined by the lifetime of the state |n), and for múltiple decay pathways, the relative 
probability for the transition |n) — > \k). Phase decoherence and population relaxation lead to 
a dissipation superoperator (represented by an iV 2 x iV 2 matrix) whose non-zero elements are 

(l-D)kn,kn = ~^ kn k ^ 71 

^D)nn,kk = +lnk k^U (8) 

^-'D)nn,nn Sn^fc 7fa' 

Population decay and dephasing allow us to overcome kinematical constraints such 
as unitary evolution to create statistical ensembles from pure states, and pure states from 
statistical ensembles, which is important for many applications such as optical pumping. 
However, there are instances when this is not desirable such as in quantum computing, 
where these effects destroy quantum information. Thus, there are situations when we would 
like to prevent decay and dephasing. A cursory glance at the quantum Liouville equation 
for coherently driven, dissipative systems @ suggests that it might be possible to prevent 
population and phase relaxation by applying suitable control fields such that 

M 

£ fUt)[H m ,p(t)}+ihL D [p(t)} = 0. (9) 

m=l 

Unfortunately, however, a more careful analysis reveals that this is not possible, in general, as 
we shall now show explicitly for a two-level system, or qubit in quantum computing parlance. 
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3. Dynamics of a 2-level system subject to control, decay and dephasing 

The Hamiltonian for a driven two-level system with energy levels E\ < E 2 is 

H[f(t)] =H + f 1 (t)H l + f 2 (t)H 2 (10) 

where H is the internal Hamiltonian and H\ and H 2 represent interaction Hamiltonians with 
independent (real-valued) control fields f\(t) and f 2 (t), 

-i " 

1 ' 

di, d 2 are the (real-valued) dipole moments for the transition and cu = (E 2 — Ei)/h is the 
transition frequency. 

We can re-write the Liouville equation in matrix form in a higher dimensional space, 
often referred to as Liouville space. Straightforward computation shows that 

j t \p(t))) = L\p(t))) = [(l/m)(L + f x {t)U + / 2 (Í)L 2 ) +í D ]\p(t))) (11) 

where \p(t))) = (pn(t), p 12 (t), p 21 (t), p 22 (t)) T 
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7 12 is the rate of population relaxation from \2) to 1 1), 721 is the rate of population relaxation 
from |1) to \2) (usually zero), and T is the dephasing rate. 

Notice that the matrix elements of the Liouville operators Li and L 2 are zero where the 
matrix elements of the dissipation operator L D are non-zero, and vice versa. Thus, no matter 
how we choose the control fields, we cannot cancel the effect of the dissipative temis. The 
best we can do is to use coherent control to implement quantum error correction schemes 
to restore decay ed/dephased quantum states to their correct vàlues. An early contribution to 
the theory of continuous feedback for such systems is contained in [p, while a more recent 
scheme of continuous quantum error correction involving weak measurements and feedback 
has been proposed in [fZp. 



4. Dynamical Semi-group and the Bloch Ball 

The more usual real vector form for p is p B = {px,2 + P2,i, i(pi,2 - P2,i), Px,\ - P2,2, px,\ + ^2,2) • 
In the general case this is referred to as the coherence vector^. Since for the systems under 
consideration we shall take p lf i + p 2 ,2 as constant (no population loss) only the first three 
components of transform under the dynamics. For pure states the norm is constant - 
thus generating motion on the surface of a 2-sphere, the Bloch sphere. For our more general 
scenario, motion takes place in the interior of this sphere, the Bloch ball. 

The family of control hamiltonians generate the Lie àlgebra u(2), in this case a 
completely controllable system[|]]. The Lie àlgebra generated by the matrices corresponding 
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to Eq.(0) acting on p B is the inhomogeneous àlgebra gl(3, R) © -R 3 . The evolution of the 
system in time is determined by exp(Lí). Noting that the eigenvalues of Í, D are (always) 
negative, the demand that our set of operators remain bounded gives in effect a semi-group, 
with only unlimited positive vàlues of t permitted. These considerations may be generalized 
to any dimensions [|]] . 

5. Conclusions 

In this note we showed how the effects of control and dissipation can be treated geometrically, 
by the movement of the coherence vector inside the Bloch ball. We incidentally noted 
that dissipative effects could not be compensated by control dynamics alone. Without 
measurements and feedback, dissipation usually forces the system into an equilibrium state. 
One can see that without any control this state corresponds to a point on the z-axis of the 
Bloch sphere; with constant controls one can show that the state converges to a point on an 
ellipse inside the Bloch ball. The coordinates of these attractors can easily be expressed as a 
function of the dissipative terms. 

Intuitively, error correction is not possible with control fields alone because the 
dissipative terms tend to pull us inside the Bloch ball and, without quantum measurement 
feedback, we can't get back to the surface of the ball because the control fields can only 
perform rotations. 

The analysis can be treated by traditional Lie algebraic methods, by use of a dynamical 
Lie àlgebra. In the specific case of a qubit treated here, this àlgebra is the inhomogeneous 
semi-direct sum gl(3,R) © R 3 . The control fields generate a rotation àlgebra within this 
larger àlgebra. Imposing boundedness conditions on the evolution of the dynamics obtained 
by exponentiation of this àlgebra leads to a semi-group description of the evolution. 

Although due to restrictions of space we have treated explicitly only the case of a two- 
level system, the methods can be generalised without difficulty to any finite level system. 
For an iV-level system, the coherence vector is essentially an N 2 — 1 real component vector 
(for population-preserving dynamics) whose motion is restricted to the iV 2 — 2-dimensional 
Bloch sphere for non-dissipative systems, and to the interior of the corresponding Bloch 
ball in general. The associated inhomogeneous real àlgebra is given by semi-direct sum 
gl(N 2 - 1, R) © i?^ 2 " 1 ), and this gives rise by means of exponentiation to the corresponding 
semi-group which determines the dynamics. Finally, for quasi-spin systems, defined by 
symmetric population decay rates, the inhomogeneous (translation-like) terms disappear. 
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